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Abstract. This note contains sufficient conditions for the probability density function of an arbitrary 
continuous univariate distribution, supported on (0, oo), such that the corresponding Mills ratio to be 
reciprocally convex (concave). To illustrate the applications of the main results, the reciprocal convexity 
(concavity) of Mills ratio of the gamma distribution is discussed in details. 

1. Introduction 

By definition (see [7]) a function / : (0, oo) — > K is said to be (strictly) reciprocally convex if x M> f{x) 
is (strictly) concave and x f(l/x) is (strictly) convex on (0, oo). Merkle [7] showed that / is reciprocally 
convex if and only if for all x, y > we have 

fll x ( 2xy\ < f_{x) + f(y) <f (x + y\ < xfjx) + yfjy) 

[ ' ' J \x + yJ- 2 - J \ 2 ) ~ x + y 

We note here that in fact the third inequality follows from the fact that the function x H> f(l/x) is convex 
on (0,oo) if and only if x i— > xf(x) is convex on (0,oo). In what follows, similarly as in [7], a function 
g : (0, oo) — > R is said to be (strictly) reciprocally concave if and only if —g is (strictly) reciprocally 
convex, i.e. if x i-> g(x) is (strictly) convex and x i-> f/(l/x) is (strictly) concave on (0, oo). Observe that 
if / is differentiable, then x i-> f(l/x) is (strictly) convex (concave) on (0,oo) if and only if x i-> x 2 fix) 
is (strictly) increasing (decreasing) on (0, oo). 

As it was shown by Merkle [7], reciprocally convex functions defined on (0, oo) have a number of 
interesting properties: they are increasing on (0,oo) or have a constant value on (0, oo), they have a 
continuous derivative on (0, oo) and they generate a sequence of quasi-arithmetic means, with the first 
one between harmonic and arithmetic mean and others above the arithmetic mean. Some examples of 
reciprocally convex functions related to the Euler gamma function were given in [7] . 

By definition (see [9]) a function / : (0, oo) — > R is said to be completely monotonic, if / has derivatives 
of all orders and satisfies 

(-l) w / (n) (aO>0 

for all x > and n £ {0, 1, . . . }. Note that strict inequality always holds above unless / is constant. It is 
known (Bernstein's Theorem) that / is completely monotonic if and only if [9, p. 161] 

/>oo 

f{x) = / e~ xt duit), 
Jo 

where v is a nonnegative measure on [0, oo) such that the integral converges for all x > 0. An important 
subclass of completely monotonic functions consists of the Stieltjes transforms defined as the class of 
functions g : (0, oo) — >• R of the form 

dv(t) 



gi x ) 



x + t 



where a > and v is a nonnegative measure on [0, oo) such that the integral converges for all x > 0. 

It was pointed out in [7] that if a function h : (0, oo) — > R is a Stieltjes transform, then —h is reciprocally 
convex, i.e. h is reciprocally concave. We note that some known reciprocally concave functions comes 
from probability theory. For example, the Mills ratio of the standard normal distribution is a reciprocally 
concave function. For this let us see some basics. The probability density function ip : R — > (0, oo), the 
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cumulative distribution function $ : R — > (0, 1) and the reliability function $ : M — > (0, 1) of the standard 
normal law, are defined by 

1 



ip(x) = 



-x 2 /2 



/2tt 

$(x) = / <p(t)dt 

J —OO 

and 

POO 

®(x) = 1 - ®(x) = / p(i)dt. 

J X 

The function to : R — > (0, oo), defined by 

¥(x) 



to(x) 



/>oo />oo 

/ e- t2 / 2 dt = / e- xt e~ t2 / 2 dt, 

Jx J0 



x 2 /2 

ip(x) 

is known in literature as Mills' ratio [8, Sect. 2.26] of the standard normal law, while its reciprocal 
r = 1/to, defined by r(x) = l/m(x) = ip(x)/&(x), is the so-called failure (hazard) rate. For Mills' ratio 
of other distributions, like gamma distribution, we refer to [6] and to the references therein. 

It is well-known that Mills' ratio of the standard normal distribution is convex and strictly decreasing on 
R, at the origin takes on the value m(0) = y/w/2. Moreover, it can be shown (see [2]) that x m'(x) /m{x) 
is strictly increasing and x n- x 2 m'(x) is strictly decreasing on (0, oo). With other words, the Mills ratio 
of the standard normal law is strictly reciprocally concave on (0, oo). Some other monotonicity properties 
and interesting functional inequalities involving the Mills ratio of the standard normal distribution can 
be found in [2]. The following complements the above mentioned results. 

Theorem 1. Let m be the Mills ratio of the standard normal law. Then the function x H- m(\fx)/ \fx is 
a Stieltjes transform and consequently it is strictly completely monotonic and strictly reciprocally concave 
on (0, oo). In particular, if x,y > 0, then the following chain of inequalities holds 




In each of the above inequalities equality holds if and only if x = y. 

Proof. For x > the Mills of the standard normal distribution can be represented as [5, p. 145] 

/oo />oo 

From this we obtain that 



r 

Jo 



-ds, 



V x V 27T Jo x + s 

which shows that the function x i-> m{ypx)l ^/x is in fact a Stieltjes transform and owing to Merkle [7, p. 
217] this implies that the function x ^ —m(y/x)/y/x is reciprocally convex on (0,oo), i.e. the function 
x i y m(y/x) I \fx is reciprocally concave on (0,oo). The rest of the proof follows easily from (1.1). We 
note that the strictly complete monotonicity of the function x ^ m(y/x)/y/x can be proved also by using 
the properties of completely monotonic functions. Mills ratio m of the standard normal distribution is 
in fact a Laplace transform and consequently it is strictly completely monotonic (see [2]). On the other 
hand, it is known (see [9]) that if u is strictly completely monotonic and v is nonnegative with a strictly 
completely monotone derivative, then the composite function mod is also strictly completely monotonic. 
Now, since the function to is strictly completely monotonic on (0, oo) and x i-> 2(y/x)' — l/\/x is strictly 
completely monotonic on (0, oo), wc obtain that x i-> m{y/x) is also strictly completely monotonic on 
(0,oo). Finally, by using the fact that the product of completely monotonic functions is also completely 
monotonic, the function x M> m{y/x)/y/x is indeed strictly completely monotonic on (0, oo). □ 

Now, since the Mills ratio of the standard normal distribution is reciprocally concave a natural question 
which arises here is the following: under which conditions does the Mills ratio of an arbitrary continuous 
univariate distribution, having support (0,oo), will be reciprocally convex (concave)? The goal of this 
paper is to find some sufficient conditions for the probability density function of an arbitrary continuous 
univariate distribution, supported on the semi-infinite interval (0, oo), such that the corresponding Mills 



MILLS' RATIO: RECIPROCAL CONVEXITY AND FUNCTIONAL INEQUALITIES 



3 



ratio to be reciprocally convex (concave). The main result of this paper, namely Theorem 2 in section 
2, is based on some recent results of the author [3] and complement naturally the results from [2, 3]. 
To illustrate the application of the main result, the Mills ratio of the gamma distribution is discussed in 
details in section 3. 

We note that although the reciprocal convexity (concavity) of Mills ratio is interesting in his own 
right, the convexity of the Mills ratio of continuous distributions has important applications in monopoly 
theory, especially in static pricing problems. For characterizations of the existence or uniqueness of global 
maximizcrs we refer to [4] and to the references therein. Another application can be found in [10], where 
the convexity of Mills ratio is used to show that the price is a sub-martingale. 



2. Reciprocal convexity (concavity) of Mills ratio 

In this section our aim is to find some sufficient conditions for the probability density function such 
that the corresponding Mills ratio to be reciprocally convex (concave). As in [3] the proof is based on 
the monotone form of l'Hospital's rule [1, Lemma 2.2]. 

Theorem 2. Let f : (0, oo) — > (0, oo) be a probability density function and let lo : (0, oo) — > R, defined 
by uj(x) — f'(x)/f(x), be the logarithmic derivative of f. Let also F : (0, oo) — > (0,1), defined by F(x) = 
J°° f(t)dt, be the survival function and m : (0, oo) — >■ (0, oo), defined by m(x) — F(x)/f(x), be the 
corresponding Mills ratio. Then the following assertions are true: 

(a) If f(x)/ui(x) — > as x — > oo, lo' /lo 2 is (strictly) decreasing (increasing) on (0, oo) and the function 

x 3 lj'(x) 
xuj 2 (x) — xuj'(x) — 2u)(x) 

is (strictly) increasing (decreasing) on (0, oo), then Mills ratio m is (strictly) reciprocally convex 
(concave) on (0, oo). 

(b) If xf{x)/{\ — xui(x)) — > 0, f(x)/uj(x) — »• as x — > oo, lo'/lo 2 is (strictly) decreasing (increasing) 
on (0, oo) and the function 

x 2 lo'(x) -XLo(x) +2 
x i — y — — 

xlo 2 (x) — xuj'(x) — 2uj{x) 

is (strictly) increasing (decreasing) on (0, oo), then Mills ratio m is (strictly) reciprocally convex 
(concave) on (0, oo). 

Proof, (a) By definition Mills ratio m is (strictly) reciprocally convex (concave) if m is (strictly) concave 
(convex) and x i-> m(l/x) is (strictly) convex (concave). It is known (see [3, Theorem 2]) that if f(x)/uj(x) 
tends to zero as x tends to infinity and the function lo' /lo 2 is (strictly) increasing (decreasing), then m is 
(strictly) convex (concave). Thus, we just need to find conditions for the (strict) convexity (concavity) of 
the function x n- m(l/x). This function is (strictly) convex (concave) on (0, oo) if and only if the function 
x i ^ x 2 m'(x) is (strictly) increasing (decreasing) on (0, oo). On the other hand, observe that Mills ratio 
m satisfies the differential equation 

m'(x) — —Lu(x)m(x) — 1. 

Thus, by using the monotone form of l'Hospital's rule (see [1, Lemma 2.2]) to prove that the function 

(F(x) + f(x)/u(x)) - lim (F(x) + f{x)/uj{x)) 
x i ^ x m (x) = — - 



f(x)/{x 2 Lo(x))— lim f{x)/(x 2 Lo{x)) 

x—>oo 

F(x) + f(x)/u(x) 



f(x)/(x 2 uj(x)) 

is (strictly) increasing (decreasing) on (0, oo) it is enough to show that 

(F(x) + f{x)/io{x))' _ x 3 lj'{x) 



X — - 



(f{x)/(x 2 Lo{x)))' ' xuj 2 (x)-xu'(x)-2u(x) 



is (strictly) increasing (decreasing) on (0, oo). 

(b) Observe that according to [7, Lemma 2.2] the function x h-> m(l/x) is (strictly) convex (concave) 
if and only if x xm(x) is (strictly) convex (concave) on (0, oo). Now, by using the monotone form of 
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l'Hospital's rule (see [1, Lemma 2.2]) the function 

, , ,,, , . , . , > F(x) — xf(x)/(l — xu)(x)) 

x h> (xm (x))' = m(x) -x- xu(x)m(x) = y J ) " K , ,, 

(F(x) - xf(x)/(l - xlo(x))) - Hm^ (~F(x) - xf(x)/(l - xlo(x))) 
= (/(x)/(l - xc(x))) - lim (/(x)/(l-xc(x))) 

x— 7-OC 

is (strictly) increasing (decreasing) on (0, oo) if the function 

(F(x) - xf(x)/(l - xlo(x)))' x 2 uj'(x) - xu(x) + 2 



X h-> 



(f{x)/(l-xw(x)))' " xcj 2 (x)-xuj'(x)-2u J {x) 



is (strictly) increasing (decreasing) on (0, oo). Notice that we used the fact that if xf(x)/(l — xuj(x)) — > 
as x — > oo, then f(x)/(l — xuj(x)) — > as x — > oo. □ 

We note here that the reciprocal concavity of the Mills ratio of the standard normal distribution can 
be verified easily by using part (a) or part (b) of Theorem 2. More precisely, in the case of the standard 
normal distribution we have oj(x) = —x, co'(x) = —1. Consequently <p(x)/u)(x) = —<p(x)/x — >• as 
x — > oo, the function x i-> w' (x) / 'uj 2 (x) = —1/x 2 is strictly increasing and 

x 3 uj'(x) x 2 

x i y 



xlo 2 (x) — xuj'(x) — 2u>(x) x 2 + 3 

is strictly decreasing on (0,oo). This is turn implies that by using part (a) of Theorem 2 the Mills ratio 
of the standard normal distribution is strictly reciprocally concave on (0, oo). 

Similarly, since tp(x)/(l + x 2 ) — > 0, X(p(x)/(1 + x 2 ) — > 0, —(p(x)/x — > as x — > oo, the function 
x i ^ cj'(a;)/a; 2 (x) = —1/x 2 is strictly increasing and 

x 2 uj'(x) - xuj(x) + 2 2 
xuj 2 (x) — xcj'(x) — 2cj(x) x 3 + 3x 

is strictly decreasing on (0,oo), part (b) of Theorem 2 also implies that the Mills ratio of the standard 
normal distribution is strictly reciprocally concave on (0, oo). 

Thus, Theorem 2 in fact generalizes some of the main results of [2]. 

3. Reciprocal convexity (concavity) of Mills ratio of the gamma distribution 

The gamma distribution has support (0, oo), probability density function, cumulative distribution 
function and survival function as follows 



f(x) = f(x;a) 



x a - x e 



and 



T(a) ' 

F(x) = F(x; a) = 3^ = f t^e^dt 

F(x)=F(x;a) = = f°° t^e^dt, 

T(a) T{a) J x 



where T is the Euler gamma function, 7(-, •) and T(-, •) denote the lower and upper incomplete gamma 
functions, and a > is the shape parameter. As we can see below, the Mills ratio of the gamma 
distribution m : (0, oo) — > (0, oo), defined by 

/ \ / \ r(o;, x) 

mix) — mix: a) = : , 

x Q— e~ x 

is reciprocally convex on (0, oo) for all < a < 1 and reciprocally concave on (0, oo) for all 1 < a < 2. 
In [3] it was proved that if a > 1, then the Mills ratio m is decreasing and log-convex, and consequently 
convex on (0, oo). We note that the convexity of Mills ratio of the gamma distribution actually can be 
verified directly (see [10]), since 

a-l 



m(x) 



poo / f\ Q ~ 1 r°° 
J l^-J e x - t dt = J xu^e^-^du 
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and 



m'(x) = J™ ({a - l)u a - 2 ) ((1 - u)e {1 - u ^ du 

/OO pOO 
M a-l e (l-n)x du + J ^-1(1 _ u ) e (l-«)x dU) 



where the last equality follows from integration by parts. From this we clearly have that 

p OO 

m"(x)= / {a-l)(l-u) 2 u a - 2 e^- u)x du 



and consequently m is convex on (0, oo) if a > 1 and is concave on (0, oo) if < a < 1. The concavity of the 
function m can be verified also by using [3, Theorem 2]. Namely, if let u>(x) = f'(x) / f(x) = (a—l)/x — l, 
then f(x)/uj(x) tends to zero as x tends to infinity and the function x i-> cj'(a;)/a; 2 (a;) = (1 — a) /{a — 1— x) 2 
is decreasing on (0, oo) for all < a < 1. Consequently in view of [3, Theorem 2] m is indeed concave on 
(0,oo) for all < a < 1. 

Now let us focus on the reciprocal convexity (concavity) of the Mills ratio of gamma distribution. Since 

x 3 oj'(x) (1 - a)x 2 



we obtain that 



xoj 2 (x) — xlj'(x) — 2u>(x) (a — 1 — x) 2 + 2x + 1 — a' 

xV(x) V _2(a — l)(a — 2)(x 2 + (1 — a)x) 



k xuj 2 (x) — xu>'(x) — 2oj(x) J ((a — 1 — x) 2 + 2x + 1 — a) 2 
This last expression is clearly positive if < a < 1 and x > 0, and thus, by using part (a) of Theorem 2 
we conclude that Mills ratio m is reciprocally convex on (0, oo) for all < a < 1. 
Similarly, since 

arV(x) - xw(x) + 2 .x 2 + 2(2 - a) a: 

a;w 2 (x) - xw'(x) - 2uj(x) ~ x 2 + 2(2 - a)x + (a - l)(a - 2) ' 

we get 

x 2 w'(a;) +2 V 2(a - l)(a - 2)(x + 2 - a) 



^cj 2 (x) - arw'(ar) - 2w(o;) / (a; 2 + 2(2 - a)a; + (a - l)(a - 2)) 2 

and this is negative if 1 < a < 2 and x > 0. Consequently, by using part (b) of Theorem 2 we get that 
the Mills ratio of the gamma distribution is indeed reciprocally concave for 1 < a < 2. Here we used that 
if x tends to oo, then the expressions f(x)/co(x) and xf(x)/(l — xui(x)) tend to 0. 

Finally, we note that the convexity (concavity) of x m(l/x) can be verified also by using the integral 
representation of Mills ratio of the gamma distribution. More precisely, if we rewrite m(x) as 



m(x) 

then 



^( 1 + H)"-e-d„ 
x 2 m'(x) = -J ( a -1)^1 + -J ue- u du 



and 

r°° / n \ "-3 7/ 2 

[x 2 m'(x)]' = / (a-l)(a-2)(l + -) -^e~ u du. 
Jo v x > x 

This shows that x n- x 2 m'{x) is decreasing on (0, oo) if 1 < a < 2 and increasing on (0, oo) if < a < 1 

or a > 2. Summarizing, the Mills ratio of the gamma distribution is reciprocally convex on (0, oo) if 

< a < 1 and reciprocally concave on (0, oo) if 1 < a < 2. When a > 2 £/ie functions x t— > m(x) and 

x i y m(l/x) are convex on (0, oo), thus in this case m is nor reciprocally convex and neither reciprocally 

concave on its support. 
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